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Abstract 

We study the behaviour of the electronic chemical potential in YBa2Cu306+5 
near the superconducting transition using the spin-wave exchange theory of 
pairing. We find that the experimental value of the jump in the tempera- 
ture derivative of the chemical potential is inconsistent with the theoretical 
value calculated under the assumption of constant hole concentration. It is 
suggested that charge transfer between chains and planes takes place and 
corresponding calculations are done. 
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The temperature dependence of the chemical potential near the superconducting transi- 
tion has attracted some attention recently It has been shown in several models that 
there is a jump in the temperature derivative of the chemical potential. In Ref. in a 
simple BCS model with quadratic dispersion and s-wave pairing, the relationship between 
the chemical potential in the superconducting state /Xg and the chemical potential in the 
normal state fin was found to be 

Hs = f^n- -7—, (1) 



where A is the gap at the Fermi surface. Near the critical point A oc ^JT^ — T, and from 
Eq. (|l]) it follows that the temperature dependence of the chemical potential is not smooth. 
This effect has been observed experimentally [§. Measurements of the work function of the 
YBa2Cu306+5 compound with the critical temperature of 90 K gave 



C^(/in - /is 



dT 



= 0.12 ±0.02. (2) 



(Throughout the paper, the Boltzmann constant ks = 1.) 

The purpose of the paper is to present a quantitative study of this effect specifically 
for YBa2Cu306+5. Using the magnetic fluctuation exchange theory of high temperature 
superconductivity, we find a value of A{dfi^/dT) much larger than the experimental value @ 
and we show that the discrepancy may be explained by taking into account the presence of 
chains. 

Before calculations, let us discuss various contributions into the temperature dependence 
of the chemical potential. The quantity measured in the experiment |Q is the single electron 
work function. To analyse it we must consider the single particle Green's function. The 
mean field approximation corresponds to the account of the leading order anomalous self- 
energy (Fig. 1) in the Dyson equations for the normal and anomalous Green's functions [§. 
With a normalisation condition, it gives an effect as in Eq. (|I]). In the next order, there is the 
self-energy contribution presented in Fig. 2, which is due to modification of the spin- wave 
(wavy line) at the transition point. Analysis of this diagram is very complicated. At zero 
temperature, this problem was recently studied in [|l^]. We plan to analyse the problem 



more generally in the future. In the present work we neglect the contribution of the diagram 
in Fig. 2 i.e. we stay within the mean field approximation. 

Another aspect to consider is the external influence on the electronic subsystem. Equa- 
tion (|1|) corresponds to the usual BCS mean field theory with constant electronic density. 
The experiment corresponds to constant external pressure, however due to interaction with 
the lattice, the electronic subsystem is neither at constant density nor at constant pressure. 
[It appropriate to note here that, at constant electronic pressure, according to Landau the- 
ory of second order phase transitions there is no discontinuity in dfi/dT (see discussion in 
Ref. 0).] We thus have to understand first which theoretical value should be compared 
with the experiment. 

There is a rearrangement of the lattice under the superconducting transition. Its influence 
on the single electron Green's function is described by the diagram in Fig. 3, with the dash 
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line representing phonons. Since there are two anomalous Green's functions in the loop 
this contribution is proportional to . However, it can be neglected if the interaction with 
phonons is weak and the lattice is hard, We therefore take the electron desity to be constant. 

Let us now outline the mechanism of pairing induced by spin- wave exchange in the t-J 
model |TT],IT2[. The t-J model is defined by the Hamiltonian 



^ = E (4.rf„.. + H.c.) + J 5: S„ ■ S™, (3) 



<nm>(T <n'm> 



where dl^^ is the creation operator of a hole with spin cr (cr =t, i) at site n of a two- 
dimensional square lattice and < nm > are nearest-neighbour sites. The dj^^ operators act in 
the Hilbert space with no double electron occupancy. The spin operator is S„ = \d'l^cra/3dni3. 
Below we set J and the lattice spacing equal to unity. 

At half- filling (one hole per site), the t-J model is equivalent to the Heisenberg anti- 
ferromagnetic model which has long-range order in the ground state. We denote the wave 
function of this (undoped) ground state by |0) and use it as the starting point for considering 
the doped state. It is well known that doping destroys long-range antiferromagnetic order. 
However, we may use |0) and the corresponding set of excitations as the basis in the problem 



with doping. The effective Hamiltonian for the t-J model as derived in [p!3HT5|| is 



HcS = J2 ^k^L^ka + J2 ^q("i"q + /^J/^q) + ^h.sw + ^^h,h- (4) 



It is expressed in terms of the usual spin waves on antiferromagnetic background and /3q 
and composite hole operators /iko- {cr = ±1/2 is the pseudo-spin index). The summations 
over k and q are over the Brillouin zone of one sublattice (defined by |A;^| + < tt). 

The wave function of a composite hole can be represented as ^/'ko- = ^ko-|0). At large t 
the operator /ij^^ has a complicated structure. At t = 3, the weight of the bare hole in ^/^ko- 
is about 25%, the weight of configurations "bare hole + 1 magnon" is ~50%, and the rest 
are more complicated configurations. The dressed hole is a normal fermion. For t ^ 5 the 



dispersion can be approximated by the expression 



Ck = const -^0.43 + 4.6 12- 2.8 t^^a + i^2(cos A:^ - cos fcj,)^ (5) 

where 7k = ^{coskx + cos ky). We use the rigid band approach, i.e. we always use the 
dispersion Eq. (^ which was derived for the undoped model. The basis for this approach 
is that the results are not highly sensitive to the dispersion and there is some evidence, 
including experimental, that the changes in the dispersion caused by doping are moderate. 
Near the band minima kg = (±7r/2, ±7r/2), the dispersion (^) can be presented in the usual 
quadratic form 

ek = ^PipI + IP2pI (6) 

where pi (^2) is the projection of k — ko on the direction orthogonal (parallel) to the face 
of the magnetic Brillouin zone, and (3i ~ 0.65 1 at t ^ 0.3. The coefficient (32 is small 
{j32 <C Pi) and the dispersion is almost degenerate along the edge of the magnetic Brillouin 
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zone (7k = 0). According to Refs. |T7|-|T^, ~ O.lt at t ^ 0.3. For the physical value of 
t 3 it means ~ 0.3. Angle resolved photoemission experiments |20] indicate a smaller 
value of — 0.1, which can be modelled by adding a small diagonal hopping term to the 
original t-J Hamiltonian. While being crucial in some experiments (such as the Hall effect), 
when pairing is considered, the exact value of (32 is not very important as long as it is small. 
In the present paper we use = and [32 = 0.3. 
The spin-wave dispersion in (H) is 



2Jl 



7q- 



(7) 



A certain approximation is made here since in the presence of mobile holes long range 
antiferromagnetic order is destroyed and a gap is formed in the spin-wave spectrum. The 
used approximation assumes that magnetic order is preserved at the scales most relevant to 
pairing. We checked that the results are not dramatically affected if a small gap is introduced 
in Eq. (0). 

The interaction of composite holes with spin waves is of the form |I^p!7|, 



h,sw~ 



^k,q {hi 



-qi^kT'^q + ^k+qT^ki/5q + H.C. 



k,q 



1/2 



(7kf/q + 7k+qVq 



(8) 



— i are the 



where N is the total number of sites, f/q = -|- | and = — sign(7q 
parameters of the Bogolyubov transformation diagonalising the spin-wave Hamiltonian. At 
t = 3, the coupling constant / in Eq. (|^) is / = 1.8 [Q. It is strongly renormalised (by a 
factor of order J/t) compared to the coupling constant for bare holes but, importantly, the 
kinematic structure of the vertex remains the same (see [p!^ , p!^ , p!8| , pl|] ) . 

Finally, there is a contact hole- hole interaction ifh,h in the effective Hamiltonian (^. 
It is proportional to some function A{t). For small t this function approaches the value 
A{0) = — |, which reflects the well known hole-hole attraction induced by reduction of the 



number of missing antiferromagnetic links. However A{t) vanishes a.tt ^ 3 |jT3|. We therefore 
neglect the contact interaction ifh,h and consider only the hole-spin-wave interaction -ffh,sw 
Because of the nearly fiat hole dispersion and extremely low Fermi energy (ei? ~ 100 K 
at optimal doping), the spin wave velocity is much higher than the typical hole velocity. It 
means that the spin-wave exchange interaction is essentially instantaneous and retardation 
can be ignored. The wave function of the superconducting state is constructed in the usual 
BCS form involving dressed holes 



i^)=nK+t^k4T^iki)io), 



(9) 



where ul. + v^. 
equation 



1 and |0) is the undoped state. For the gap we then have the usual BCS 



/^k = 

i?k = + 



Vkk' tanh — — , 



2EJ 



(10) 



ek - /i, 
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TABLE I. Partial derivatives of hole concentration at the transition point. The hole concen- 
tration is X = 0.15. 



/32 






T 


{dx/dn)n 


{dx/dfi)s 


(ax/9r)„ 


(dx/dT), 


0.0 


0.006 


0.010 


0.042 


1.97 


2.45 


0.92 


-0.55 


0.3 


0.134 


0.139 


0.039 


1.36 


1.69 


0.84 


-0.44 



with the chemical potential /i fixed by the hole concentration 

a; = y fl-iHtanh^V (H) 

The interaction in Eq. ( [TI1| ) is 

r) fl'-k,qfi'-k',q 

^ ~ , iTi — \FT ^ 

-^^q - l^kl - l^k'l, 

where q = k + k'. The diagrams representing the exchange of one spin wave are shown in 
Fig. 4. The kinematic structure of the interaction is better seen if we use Eq. (|^) and write 

4/' 27k7k' - (7k + 7g07q 

^-k,q^-k',q - ^ (1 _ ^2)1/2 " ^^'^) 

The structure of the interaction (|T2|) resulting from the hole-spin-wave vertex of Eq. ([T3|) is an 
important difference of the used approach from the phenomenological theories of supercon- 



ductivity ||22| which take the hole-spin-wave vertex as a constant without any k-dependence. 
We believe that Eq. (|T2|) is more appropriate in the context of strongly correlated systems 
since it is derived from the microscopic Hamiltonian. 

Details of an approximate analytical and exact numerical solutions of Eq. (^) are dis- 
cussed in [|ri|] and [|12|. There is an infinite set of solutions with different symmetries but 



the strongest pairing is in the dx2-y2-wave channel. 

The calculated chemical potential for hole concentration x = 0.15 is presented in Fig. 5. 
To study analytically the temperature dependence of the chemical potential, it is convenient 
to use partial derivatives of concentration. For the jump in the temperature derivative of 
the chemical potential we have 



d{l2n - /is 



dT 



T=T, 



{dx/dT)s _ jdx/dT)^ 
{dx/dfi)s {dx/d^)n ' 



where indices s and n refer to derivatives being taken below and above the transition point, 
respectively. In experiments, electron rather than hole chemical potential is measured and 
the appropriate quantity is therefore A^dfi'^/dT) = —A{d^/dT). 
The results of our calculations are presented in Table I. 

Substituting data from Table I into Eq. (|T^) or from Fig. 5, we find A{dfj,ydT) = 0.69 
for (32 = and A{dfi'^/dT) = 0.88 for = 0.3, which is more than 5 times larger than the 
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experimental value (0). In the rest of the paper we examine how the presence of chains may 
affect this result. 



The chains in YBa2Cu306+<5 are known to be charge reservoirs and we can expect that the 
presence of a reservoir will reduce variations of the chemical potential (in the extreme case 
of an infinite reservoir the chemical potential is always constant). The conserved quantity 
is the total hole concentration 5 = 2x + x', where x is the in-plane hole concentration, x' is 
the hole concentration in the chains, and we have taken in account that there is one chain 
for two planes. The expression for the jump in dfi/dT is now 



dT 



T=Tc 



2{dx/dT)s + jdx'/dT) 2{dx/dT)n + jdx'/dT) 
2{dx/di2)s + (dx'/dfi) ~ 2(9x/9/i)„ + {dx'/dfi) ' 

The number of holes in the chain is found as in the usual filling of a band, 

dk 



(15) 



(16) 



where n^(e) = exp(^^) + 1 and e'^ is the dispersion of the chain band. Below the 
transition, there may be superconducting order in the chains induced by proximity effect. 
However, the gap is expected to be small, especially in the vicinity of the transition point 
and in any case it would have no effect on the dependence of x' on /i. 
We will take the chain dispersion in the simple quadratic form 



el = eo + kPk^ 



(17) 



where eo is the chain band minimum relative to the bottom of the plane band and (3 is the 
inverse effective mass. Then at low temperatures (T ^ fi — cq) 



X 



23/2 



TT 



(/^ - eo) 



24 (/X - eo) 



+ 



(18) 



It is important to note here that the shift eo between the plane and chain bands is determined 
by self-consistent distribution of Coulomb charge and is a function of x'. Differentiating ( p!8[ ) 
we find 



and 



' dx'' 



' dx' 



1 

TT 



1 



(19a) 



'dx^ 
dT) 



1 + /.^ 



deo 
' dx' 



-1 



TT 

6 



T 



16T 



From Eqs. (19), we have the relationship 



4T 

'^^{x'f 



(3{n- eo)3/2 37r2/32(a; 



' dx'^ 



'\3 



(19b) 



(20) 
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which allows us to eliminate the unknown dependence eo{x'). To produce the experimental 
value of A{d^'^/dT), we need to take {dx'/dfi)T = 14. For any reasonable value of f3, the 
contribution of (dx'/dT) from Eq. ( |19b| ) into Eq. ( ]I5| ) is negligible (we take x' = 0.7). 

Now we are in a position to estimate the amount of hole transfer between planes and 
chains. Using {dx'/dfi)T = 14 and Eq. (^), we have 

Ax' = x'(T,) - x\0) = 14 (^A/i - ^1^^ , (21) 

where A/i = /i(Tc) — /i(0). For the in-plane concentration, we have 

Ax = 2.45A/i - 0.27 (22a) 

for /?2 = and 

Ax = 1.69A/i - 0.22 (22b) 
for = 0.3. Now using the condition 2 Ax + Ax' = 0, we find 

Ax = 1.21— ^-0.20 T, (23a) 

for (32 = 0.0 and 

Ax = 0.91^- 0.18 (23b) 

for (32 = 0.3. For further estimates we will take the average of the two expressions. The 
value of Ax depends strongly on (3. Ax is positive (holes leave the chains with increasing 
temperature) if (3 < Per = 0.45. In real units (using J = 0.15 eV and the lattice constant 
a = 3.8 A), it corresponds to effective mass m* > m*j. = 7.5me. As an example, for (3 = 0.20 
(m* = 17me) the change in the in-plane concentration is Ax = 0.009, so that the in-plane 
hole concentration at is larger by 6 % compared to T = 0. 

To summarise, we have studied the behaviour of the chemical potential near the super- 
conducting phase transition using the magnetic fluctuation exchange theory of superconduc- 
tivity. We suggested that the found discrepancy with the experimental value of the jump 
in the temperature derivative of the chemical potential is be caused by changes in the in- 
plane hole concentration. We presented a model which shows that even a small hole transfer 
between chains and planes significantly reduces the value of the jump. Further experiment 
data (involving other superconductors) would lead to better understanding of the problem. 

We are grateful to D. van der Marel for drawing attention to the problem and to D. I. 
Khomskii and M. Yu. Kuchiev for stimulating discussions. 
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FIGURES 

FIG. 1. The diagram for the leading order anomalous self-energy. 

FIG. 2. The influence of magnon self-energy on the single electron Green's function. 

FIG. 3. The influence of the lattice on the single electron Green's function. 

FIG. 4. The diagrams for exchange of one spin wave. 

FIG. 5. The chemical potential //. as a function of temperature T. The dash line is for 
P2 — and the sohd line is for P2 — 0.3. The hole concentration is constant, x — 0.15. 
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Fig. 4 



FIGURE 5 
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